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ABSTRACT 

We develop an empirical model to estimate mass-loss rates via coronal mass ejections (CMEs) for 
solar-type pre-main-sequence (PMS) stars. Our method estimates the CME mass-loss rate from the 
observed energies of PMS X-ray flares, using our empirically determined relationship between solar 
X-ray flare energy and CME mass: log(McME [g]) = 0.63 x log(£ , nare [erg]) — 2.57. Using masses 
determined for the largest flaring magnetic structures observed on PMS stars, we suggest that this 
solar-calibrated relationship may hold over 10 orders of magnitude in flare energy and 7 orders of 
magnitude in CME mass. The total CME mass-loss rate we calculate for typical solar-type PMS stars 
is in the range 10 -12 -10 -9 M Q yr _1 . We then use these CME mass-loss rate estimates to infer the 
attendant angular momentum loss leading up to the main sequence. Assuming the CME outflow rate 
for a typical ^1 M© T Tauri star is < 10~ 10 Mg yr -1 , the resulting spin-down torque is too small 
during the first ~1 Myr to counteract the stellar spin- up due to contraction and accretion. However, 
if the CME mass-loss rate is > 10 -10 M© yr _1 , as permitted by our calculations, the CME spin-down 
torque may influence the stellar spin evolution after an age of a few Myr. 

Subject headings: Stars:activity — Stars:evolution — Stars: mass-loss — Stars: pre-main sequence — 
Stars: solar-type 
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1. INTRODUCTION 

Solar-type pre-main-sequence (PMS) stars typically 
evince magnetic activity in a variety of forms, including 
strong and ti me- variable coronal emission at X-r ay wave- 
lengths (e.g. lGudelll200l iPreibisch et al.ll2005l ). strong 
and time-variable c hromospheric emission at ultravio- 
let wavelengths (e.g . lYang et al.ll2012t IPetrov et al]|2011t 
iValenti et al.|[l993|) and in chromospheric emission lines 
(such as Ha; e.g. iWhite fc~B asri 20031 iFernandez et al.l 
ll995UWalter fc Kuhilll981l). large spot-modulated photo- 
metr i c variability (e.g. | Herbst et "al"l l2007; St assun et al.l 
IT999llHerbst et aLHltiM T and non-thermal polarized ra- 
dio continuum emission is a classic indicator of mag- 
netic activity (e.g., the PMS star Hubble 4, iSkinnerl 
1993). These observations have motivated direct Zee- 
man measure ments revealing fiel d strengths of typically 
a few kG (e.g. I Johns-K rulll 120071 iR einers fc Basrt liM 
IJohns-Krull et al.l I2004D . The strong magnetic fields as- 
sociated with this ubiquitous activity are thought to be 
central to a number of key physical processes in PMS 
stars, including specifically the transfer of mass and of 
angular momentum from and to the circumstellar envi- 
ronment. 

X-ray activity observed in PMS stars bears a striking 
resemblance to solar X-ray activity, albeit scaled up by 
several orders of m agnitude in both energy and frequency 
of occurrence le.g. JPeres et alj2004trGetrnan et aLll2005l; 



IPeres et all 12005 lAudard et alJl200l . Generally, PMS 
X-ray flares are well described by standard solar flare 
models (e.g., iReale et"al~l 119971 ). though in the case of 
some exceptionally long-lived flares it has been neces- 
sary tc^_ad<ij3elf^clip^ the star 
(i.e JSkinner et al.|[l997l : I Johnstone et alj2012j) . In addi- 
tion, with the exception of a few very strongly accreting 
systems where the X-ray variabilit y is at least in part 
accretion driven (e.g., V1647 Ori; iKastner et all I2004t 
lHamaguchi et al.l I2012D . large-scale studies of simulta- 
neous X-ray and optical variability in PMS stars show 
that the X-ray variability is p r incipa lly caused by coronal 
activity (jStassun et al.l 120061 120071 ). Young, solar- type 
stars appear, then, to have flaring magnetic field con- 
figurations that behave like, and can be interpreted as, 
scale d-up solar type flaring fields (see also lFeigelson et all 
20071 for a recent overview). 

During the early, active accret ion phase of a PMS star's 
evolution (first few Myr; e.g. Qjaisch ct al. 200lJ), the 
strong stellar field may be important in governing the 
interaction of the star with its protoplanetary disk. This 
magnetic star-disk interaction likely controls the funnel- 
ing o f circumstellar m aterial from the disk onto the star 
(e.g. iShu et al"1ll994f) . is likely central to the launching: 
of ou t flows including jets and winds (e.g. IHavashi et al.l 
119961: IMatt fc Pudritzl I2005al) . and therefore is likely 
a key mediator in the net flux of mass and angular 
momentum onto and away from the young star (e.g. 
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Matt fe Pud"rit9l2005al iMatt et alj|2010t lOrlando et all 
2011). 

Even after the phase of active accretion subsides, 
scaled-up solar-type winds and elevated solar-type X- 
ray activity do not decline t o present-day solar lev- 
els until approximately 1 Gyr ((Wood 20041 iWood et al.l 
120051: iGuedel et all 119971) . potentially affecting the on- 
going evolution of the star's angular momentum and 
the circumstellar e nvironments of young planets (e.g. 
IVidotto et aLl[20T0l ). 

Extremely large magnetic loops — extending out to 
tens of stellar radii from the stellar surface — have 
been inferred among the most powerful X-ray flaring 
PMS stars from the C handra Orion Ultra-deep Project 
(COUP) sample (e.g. jFavata et al.l 120051: iMassi et al.l 
1200% ISkellv et al.l 120081) . implying" very large magnetic 
lever arms that are potentially able to s hed s ignificant 
angular momentum (jAarnio et al.1 120091 l2012t ) . While 
it was at first appealing to picture these large mag- 
netic loops as being linked to circumstellar disks, per- 
haps either arising f rom a star-disk reconnection event 
(jOrlando et alJl20ll or sim ply maintaining loo p stabil- 
ity via anchoring to the disk. lAarnio et ail ([2010D showed 
that the majority of these stars lack massi ve disks within 
reach of the loops. Indeed, earlier work bv lGetman et al.l 
(2008) came to similar conclusions, inferring from the ap- 
parent absence of large magnetic loops among stars with 
close-in disks that perhaps it is the absence of the trun- 
cation by an inner dust disk which allows these loops to 
be so large in the first place. 

On the Sun, the eruption of such magnetic struc- 
tures can drive coronal mass ejections (CMEs), which 
are observed in Thomson -scattered optical light (cf., 
iVourlidas fe Howard! 120061 ) . While the faint scattered 
light of CMEs cannot yet be directly detected from most 
other stars, there have been r ecent observations of y oung 
stars with the FUSE satellite (Leitzi nger et al.ll2011[ ) and 
the detection of a larg e flare-associated mass ejection in 
the P MS star Z CMa (jStelzer et al.l [20091: iWhelan et all 
120101 ) . More generally, CMEs are expected from low- 
mass stars by analogy to the Sun as well as from basic 
escape velocity considerations for the flaring material; in- 
deed, CMEs are believed to have been detected on M and 
K dwarfs from balmer line asymmetries, transient spec- 
tral line absorption co mponents, and EUV dimmings (cf., 
iLeitzinger et al.ll201lL and references therein) . In partic- 
ular, the extreme X-ray flares observed on PMS stars 
are expected to h ave associated extreme CMEs (e.g., 
lAarnio et al.ll20ill) . the ramifications of which — for the 
evolution of the star, for the protoplanetary environment, 
and for the star's angular momentum evolution — depend 
on the resulting mass-loss rate and the evolution of that 
mass loss over the course of the star's PMS evolution. 

A number of studies have used emission line trac- 
ers to empirically estimate the mass-loss rates from 
PMS winds during the actively a ccreting, classical 
T Tauri star (CTTS ) phas e (see, e.g.. iKwan etaLl 120071: 
lEdwards et al.l I2006L I2003L and references therein) , and 
such accretion-powered P MS winds have been explored 
in detail theore t ically (e.g. IMatt fc Pudritzll2005aLl2 008a: 
Cranmerl |200l IMatt fc Pudritzl |2008bt iCranmerl 120091: 
Matt et al.l 12010( 7! Meanwhile, reliable empirical mea- 
sures of mass-loss rates for non-accreting PMS stars, 
particularly from impulsive CME-type events, are vir- 



tually nonexistent, and simulations of winds from non- 
accreting, weak-lined T T auri stars (WTTSs) h a ve onl y 
recently been investigated (jVidotto et al.ll2009allbl . l2010l) . 

Here we aim to make first empirical estimates of mass- 
loss rates for PMS stars resulting from CMEs, and then 
use these mass-loss rate estimates to infer the attendant 
angular momentum loss leading up to the main sequence. 
Specifically, we develop a method to estimate CME mass- 
loss rates from observed PMS X-ray flares, calibrated to 
solar observations and models, and then calculate the 
resulting braking torques in the context of PMS stellar 
evolutionary models. 

We begin by summarizing our previous work to estab- 
lish an empirical relationship between X-ray flare flux 
and CME mass for the Sun, and describe how we then 
extend that relationship to solar- type PMS stars (Sec. [2]). 
In Section[3j we use our empirical flare/CME relationship 
to derive CME mass-loss rates for PMS stars both empir- 
ically (using directly observed flare properties) and ana- 
lytically (using functional representations of flare proper- 
ties), and furthermore test our methodology on the Sun, 
verifying that we are able to recover its well measured 
CME mass-loss rate from its observed X-ray flares. Next, 
in Section |4] we calculate the attendant torque on PMS 
stars exerted by the CMEs, and assess the likelihood that 
these braking torques may serve to prevent the stars from 
spinning up as they contract toward the main sequence. 
Finally, we discuss our findings in Section [5] and summa- 
rize our conclusions in Section [6] 

2. A SOLAR-CALIBRATED RELATIONSHIP BETWEEN 
PMS X-RAY FLARES AND CME MASSES 

In lAarnio et al.l (|2011| ). we cross- matched 10 years of 
solar flare and CME observations in order to determine, 
for associated flares and CMEs, whether there are cor- 
relations between flare and CME properties. Approx- 
imately 1000 solar flares and CMEs were found to be 
spatially and temporally associated. The resulting rela- 
tionship between flare flux and CME mass is well fit by 
a single power law of the form: 

log(M [g]) = 18.7 + 0.7 x log(F [W uT 2 ]), (1) 

where M is the CME mass and F the flare flux. Flare 
flux and CME mass hold this log-linear relationship over 
a few dex in flux and mass. 

Solar flares are, by convention, classified by their peak 
X-ray flux, while for stellar X-ray flare observations, flare 
luminosities are reported. Therefore, to relate our so- 
lar flare flux/CME mass relationship to observations of 
PMS stellar flares, we convert the solar flare X-ray fluxes 
into energies, and re-frame the CME mass/flare flux re- 
lationship into a CME mass/flare energy relationship. 
To determine solar flare energy, we integrate the flux 
of the solar X-ray flare light curves from flare start to 
end time. For simplicity, we approximate the flares by 
their peak energy and assume a linear decay. Specifi- 
cally, each flare's total energy output equals one half of 
the observed peak flux, times the observed flare dura- 
tion, times 47r(lAU) 2 . Our final solar CME mass/flare 
energy relationship is shown in Figure [1] and it has the 
functional form: 

M = KmE?; (2) 
in this and all following cases, M will denote the CME 
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mass and E the flare energy. Km is (2.7±1.2)xl0 3 in 
cgs units, and /3 = 0.63±0.04. 

In adopting this relationship, we take advantage of 
the empirical correlation between associated flares and 
CMEs, but we make no assumption or statement of a 
causal relationship between hares and CMEs. Indeed, it 
is well established that not every flare is associated with 
a CME, and not ev ery CME is associated with a flare. 
lAarnio et all (j2011l) found that 13% of CMEs are flare 
associated and ~9% of flares are CME associated; this 
lack of correspondence is at least partially driven by ob- 
servational biases, and any correction for over- or under- 
estimating the number of CMEs associated with flares is 
likely less than approximately a factor of two. The most 
energetic solar flares, however, are the most likely to be 
CME-associated , and these CM Es are the most massive 
(I Andrews! [20031 lYashiro et all [20051 [2001 lAarnio et all 
I2011L to name but a few studies on the matter). Since 
young stars' flares are very powerful with respect to solar 
flares and the most energetic solar flares have the high- 
est rate of CME coincidence, we could suppose that a 
flare/CME association probability is close to unity for 
TTS flares. 

Since we cannot yet directly observe CMEs on young 
stars, we will use our solar-calibrated relationship (eq.dj 
as applied to X-ray flare observations of young stars to 
determine limits on the mass losses from young stars via 
solar-analog CMEs. However, PMS X-ray flare energies 
are 5-7 orders of magnitude stronger than solar flares, 
requiring a large extrapolation of the solar flare/CME 
relationship. Therefore, it is prudent to demonstrate that 
the likely masses of CMEs associated with the extreme 
X-ray flare energies observed on PMS stars are in fact 
consistent with the solar flare/CME relationship. To do 
this, we have estimated the masses associated with the 
32 most powerful PMS X-ray flares observed in Orion 
by COUP ([Favata et all 120051 ) . Strictly speaking, these 
masses represent the masses of post-flaring loops, which 
might not be equivalent to the CME masses, but in any 
case by analogy to the Sun should be related to the CME 
masses. 

We use the flaring loop physical parameters fo r these 
32 extreme PMS flares as derived bv lFavata et al.l (|2005[ ). 
who used the solar-calibra t ed un iform cooling loop 
(UCL) model of iReale et all ([19971 ) to infer the flaring 
loop lengths and their densities, the latter inferred from 
the flare peak emission measures. We adopt the typi- 
cally assumed cylindrical loop shape and the typically as- 
sume d loop aspect ratio(radius/length = 0.1 Rea le et all 
119971) . permitting the flaring loop volumes to be calcu- 
lated. The loop masses then follow from the volumes 
and densities. In this way, we obtai n flaring loop masses 
ranging from 10 19 g to 10 22 g for the lFavata et al.l (|2005l ) 
sample. These 32 flaring loop masses and their associated 
flare energies are shown in Figure [Tj as the point with er- 
ror bars at the far upper right, representing the mean and 
standard deviation of the sample. The large scatter of 
~1 dex in the masses is likely due in part to the simplistic 
assumptions of the UCL model upon which the mass es- 
timates are based. For example, the UCL model assumes 
that a single magnetic loop is involved in the observed X- 
ray flare, while on the Sun, arcades of s maller loops are 
often involved in a given flaring "event." [Vaananc n et al.l 
( 2009) suggest that the height of such loop arcades can 



differ from the UCL inferred loop heights by factors of 2- 
10, consistent with the scatter in loop masses we have es- 
timated in Figure [1] In any case, the extreme PMS CME 
masses so inferred are consistent to within a factor of 10 
with that predicted by the solar-calibrated flare/CME re- 
lationship. Note that we have not adjusted the fit from 
the solar data, but merely extrapolated it, which lends 
strong support to the use of this same solar-calibrated 
flare/CME relationship over the full range of solar to 
TTS flare energies. 

3. MASS LOSS VIA CMES 

Having established the form of the relationship be- 
tween flare energy and CME mass, in this section we 
derive total CME mass-loss rates for PMS stars, both 
(1) empirically using the observed distribution of flare 
energies directly (Sections [37TTT1 and l3~2~Tj) . and (2) ana- 
lytically using the observed flares to describe a functional 
form for the distribution of flare energies (Sections 13.1.21 
and !3.2.2] l. First, however, we apply our methodology to 
the solar case as a confirmation that our method recovers 
the observed CME mass-loss rate for the Sun. 

3.1. Solar case 

We start by demonstrating our method for the solar 
case as a benchmark. For the sun, we have a measured 
CME mass distribution, from which we can directly cal- 
culate a CME mass loss rate to compare to the mass-loss 
rates that we infer from the X-ray fl ares. 

In the LASCO CME database ([Gopalswamv et al.1 
[20091) . 13,862 CMEs are reported from 1996 to 2006, 
spanning almost a full solar activity period including 
minimum and maximum. Of those, 6,733 CMEs have 
well constrained mass measurements, and an additional 
1,395 have masses flagged as highly uncertain. We show 
the full rate distribution of reported CME masses (num- 
ber of CMEs per year as a function of CME mass), in- 
cluding both well-constrained and highly uncertain CME 
masses, in Figure [2] (upper middle panel). We have in- 
cluded here the additional halo CMEs (those projected 
toward the Earth, or wh ich are so wide as to appear as 
if they were) for which lAarnio et al.1 (|2011| ) estimated 
masses. 

Summing up this entire distribution gives a CME mass 
loss rate of 7.8xl0 -16 M Q yr -1 . Put another way, be- 
tween 1996 and 2006, the Sun shed at least 1.563xl0 18 
g yr -1 via CMEs. This value represents a lower limit to 
the mass loss rate via CMEs from the Sun, as just under 
half of the reported CMEs have measured masses and 
there may be C MEs that escape obs ervation, especially 
at lower masses. lAarnio et all (|201lD found 1,153 of the 
CMEs with measured masses to be associated with flares 
(the masses of these CMEs are also shown in Figure [2j 
upper middle panel) ; these CMEs were used to derive the 
relationship between flare energy and CME mass (eq. [5]) . 
Even though only ~15% of the CMEs with measured 
masses are associated with flares, these flare-associated 
CMEs constitute 40% of the total observed CME mass 
loss rate. 

The above directly measured CME mass-loss rate thus 
provides a benchmark that our flare-inferred CME mass- 
loss rate below should reproduce. To compare the CME 
mass-loss rate that we will infer from flares below to the 
observed solar case, we will find it convenient to represent 
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the solar CME mass distribution (Fig. [21 upper middle 
panel) as a differential CME mass distribution (dN/dM), 
and this is shown in Figure lower middle panel. 

3.1.1. Empirically determined solar CME mass-loss rate 

To estimate the solar CME mass-loss rate from flares, 
we require a flare energy distribution to use with our 
empirical relation between flare energy and CME mass 
(cq. [2]). In order to do this as we will below for TTSs, 
we can determine a solar flare frequency f/; r as done by 
lAlbacete Colombo et all (|2007t) for the COUP data. The 
GOES satellite recorded 22,674 solar flares from 1996 
to 2006, so the total f/; r = 0.065 ksec -1 (this value is 
0.0007 for the flares observed by the COUP; see Sec. l5"2"Tl 
below). The energy distribution of observed flare rates 
(number per year as a function of energy) is shown in 
Figure [21 top left panel. As with the CMEs, we convert 
this flare energy distribution into a differential energy 
distribution (dN/dE; Fig. [2] lower left panel) and fit a 
power law of the form: 



(3) 



our fit above an estimated completeness threshold (see 
below) of E cut =10 275 erg (vertical line in Fig. [21 lower 
left panel), gives a =1.92 (±0.02) and K E = 4.7xl0 20 
(±3%) in cgs units. Our value for a is consist ent with 
those reported in the literature: [Hudson (199lj) reports 
a = 1.8 in ^ = AW~ a , where W is the total energy 
radiated by the flare, and below the GOES detection 
limit, RHESSI (Ramaty High Energy Solar Spectroscopic 
Imager) data indicate that this distribution extends to 
lower flare energies and contin ues to exhibit powe r-law 
behavior, with a slope of ~1.5 (jChriste et al.ll2008h . 

Next, we compute the inferred mass distribution of 
CME rates by using equation [2] to convert the energy 
bins in the flare energy distribution (Fig. [21 left bottom 
panel) into mass bins, i.e., dN/dE — > dN/dM, where 

dN _ dN dE 

dM ~ ~dE X dM ( ' 

and where dE / dM is given by equation [21 The result- 
ing differential CME mass distribution is shown in Fig- 
ure [21 right bottom panel, and the corresponding non- 
differential distribution is shown in the right upper panel 
for ease of comparison with the other distributions shown 
in the top panels. The total mass in this empirically de- 
termined mass distribution is M cmc = dN/dM M = 
1.558xl0 18 g yr _1 , which impressively is equal to the to- 
tal observed LASCO solar CME mass loss rate fSec. l3.1[ ) 
to within less than a percent. 

This is reassuring that the CME mass-loss rate inferred 
directly from observed flares is able to reproduce the di- 
rectly measured CME mass-loss rate. At the same time, 
as the right column of Figure [5] shows, our inferred mass 
distribution reflects a couple of known biases: first, we 
under-predict the highest mass CMEs because our flare- 
CME relationship is not taking into account the associ- 
ation probabilities of flares and CMEs which increases 
with greater flare energy. Second, we over-predict the 
observed occurrence of lower mass CMEs likely because 
of a bias against the detection of the lowest energy flares 
(in Fig. [2j dN/dE "turns over" at lower energies). 



Thus, while our empirical flare-inferred CME mass-loss 
rate appears to work well at reproducing the directly 
measured CME mass-loss rate, we can attempt to ac- 
count for those flares and CMEs that are missed obser- 
vationally by describing the observed CME rate analyti- 
cally, as we now discuss. 

3.1.2. Analytically determined solar CME mass-loss rate 

In an effort to account for the observational biases men- 
tioned above — most importantly the missed low-energy 
flares and therefore the lowest mass CMEs — we can re- 
place the directly observed flare distribution with ana- 
lytical fits to the observationally complete parts of our 
flare energy and CME mass distributions, and assuming 
a simple power-law of the form: 
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where Mi represents our flare-inferred CME masses from 
equation [21 and 
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Equation [5] is the integral of the inferred CME mass dis- 
tribution, i.e., the integrated power law of Figure[2l bot- 
tom right panel, and Ke and /3 are from equation [2] (the 
fitted values of Km,i and 7; are reported in the plot leg- 
end). For the upper limit of the integral, we take the 
highest solar CME mass observed, M max ~6xl0 16 g. 
To determine the value of M m i n , we can either use the 
observational completeness limit of dN/dE or else sim- 
ply the minimum observed CME mass. We begin with 
the clearest observational constraint, the completeness 
limit of dN/dE. Plugging E cut (10 27 5 erg) into the flare 
energy/CME mass relationship (eq. [2]), we find a corre- 
sponding M m i n of ^10 15 g. This represents a very conser- 
vative minimum mass; clearly, the observed CME mass 
distribution contains a significant number of CMEs less 
massive than ~10 15 g. The resulting mass loss rate, then, 
as calculated with equation [5] is 6.6xl0~ 16 M yr _1 . 
Were we to simply set M m i n equal to the minimum mass 
CME inferred, ^10 13 g, we find a mass loss rate of 
2.2xl0 -15 M Q yr _1 , roughly 10% of the solar wind mass 
loss rate. 

3.2. T Tauri Star case 

In Section I57T1 we demonstrated that the observed dis- 
tribution of solar flares, when converted to CME masses 
via our empirical relation between flare energy and CME 
mass (eq. [2]) , correctly recovers the observed distribution 
of solar CMEs. We moreover developed an analytical 
representation of the solar flare distribution in an at- 
tempt to account for observational biases in the observed 
solar flare distribution, and hence in our flare- inferred so- 
lar CME mass distribution. In this section we follow the 
same procedure, now applied to PMS stellar flares in or- 
der to arrive at a PMS stellar CME mass-loss rate. 
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3.2.1. Empirically determined TTS CME mass-loss rate 

As in the solar case, we can infer an empirical 
CME mass-loss rate for TTSs by combining our solar- 
calibrated flare energy/ CME mass relationship extended 
into the TTS flare regime (eq. [5] and Fig. [IJ with an em- 
pirical flare-energy distribution function for TTS flares. 

We utilize the energy distribution of TTS flare 
rates already measured for t he COUP sample by 
lAlbacete Colombo et all (120071). Next, we convert 
the cu mulative distribution of lAlbacete Colombo et al.l 
(|2007h to dN/dE (Fig. El left panel). Proceeding as in the 
solar case (Sec. l3.1| ). this dN/dE and the CME mass/flare 
energy relationship (eq. [2] and Fig. [IJ are used to calcu- 
late a CME mass distribution (Sec. 13.1.21 Fig. [31 middle 
and right panels). Directly summing this inferred CME 
mass distribution, we obtain a CME mass-loss rate of 
Mcme = Y, dN /dMM = 6.2xl0~ 13 M Q yr" 1 . This is 
an order of magnitude more mass-loss than the present- 
day solar wind. As in the solar case, this likely represents 
a lower limit on the true TTS CME mass-loss rate due to 
observational biases against the detection of the smallest 
flares/CMEs. 

3.2.2. Analytically determined TTS CME mass-loss rates 

As in the solar case, the empirical flare-inferred Mcme 
presented above is subject to one clearly problematic 
bias: the X-ray flare detection limit. Below this limit, 
E cu _t, as in the solar case, the flare rate distribution likely 
continues to increase to ever lower flare energies; if this 
behavior is as on the Sun, the distribution likely follows a 
power law over many dex in energy. The empirical Mcme 
above, then, provides a lower limit for the mass loss. 
Here we present an analytical solution that represents 
an upper limit by attempting to account for the under- 
detected, lower energy/lower mass flare/CME events. 

The power law fit to the COUP flare energy /event 
rate distribution above the completeness limit of E cut = 
10 35 6 erg is as described in e quation [3j where a —2.1 
(jAlbacete Colombo et al.ll2007[ ) . We do not have directly 
observed bounds on TTS CME masses, so we instead per- 
form our distribution integration over flare energy and 
using our relation between flare energy and CME mass 
(eq. [2]) . The mass loss rate is then (equivalent to eq. [5l 
with substitutions): 



ME : 



K E K M 



(9) 



j\ ^ max mm J \^ ) 



(/3-a + l) 

To place a lower limit on the value of E mm (and thus 
an upper limit on Mqme), we find the energy at which 
the X-ray luminosity in flares (Lfl aro ) is equal to the total 
stellar X-ray luminosity, L x ■ 

f E ™** dN K E , , 



~ E min) = L X- 



(11) 

For the stellar Lx, we take a fiducial va lue of 2.5 x 10 30 
erg s" 1 from the empirical relationship of Prcibisc h et al.l 
(2005) for a 1 M Q star. We note that the scatter in Lx is 
rather large (see figure 3 of iPreibisch et al.ll2005l) : the la 



scatter about their linear regression fit to Lx as a func- 
tion of mass is 0.65 dex. Thus our choice of fiducial Lx 
could be as low as 5xl0 29 erg s _1 or as high as lxlO 31 
erg s _1 . We choose the maximum energy in equation 
[TT1 to simply be the most energetic flare observed in the 
COUP sample, ^10 37 erg. Evaluating equation QTJ we 
find E min > 3.4xl0 29 erg. With these parameters, we 
find an upper limit to the mass loss rate for TTS via 
equationll M CM E < 3.2xl0" 10 M© yr" 1 . 

This mass loss rate is most sensitive to the flare fre- 
quency/energy distribution slope a; were we to vary 
a within the reported uncertainty, the mass loss rates 
would change by a couple orders of magnitude: for a 
values of 2.05 and 2.2, M C me is 1.9xl0~ 9 Mg yr" 1 and 
7.1xl0~ u M Q yr -1 , respectively. This is consistent with 
the solar case and our expectation that the CME mass 
loss rates would not exceed typical wind mass loss rates 
for young stars. 



3.3. CME mass loss rates compared to steady wind, 
observed outflow rates 

We have devised a method to infer the CME mass- 
loss rate from the observed flaring rate together with our 
empirical relation between flare energy and CME mass. 
Testing this method on the Sun, our empirical method 
predicts the CME mass loss rate within less than a per- 
cent of the observed value. We then applied the same 
method using an analytical power-law representation of 
the observed flare distribution in order to attempt to ac- 
count for observational biases against the lowest-energy 
flares (and thus against the lowest mass CMEs). This an- 
alytical method predicts 2.8 times more mass loss than is 
observed. With r espect to the solar wind mass loss rate 
(10~ 14 M Q yr" 1 , lLl ll999f) . the empirical result (as well 
as the observed CME mass loss rate) is 4% of the wind 
mass loss rate, while the analytical result is that CMEs 
lose mass at 10% the solar wind rate. Given that the ob- 
served CME mass loss rate is itself underpredicted (only 
half of the CMEs reported in the database have measured 
masses), and the observational bias against the lowest 
mass flares is present in our empirical result as well as 
the actual solar data, we would anticipate the analytic 
calculation to be greater than the observed mass lost. As 
such, these two approaches — empirical and analytical — 
serve as lower and upper limits, respectively. 

In the TTS case, our empirically calculated lower 
limit on the stellar CME mass loss rate is ~6xl0~ 13 
M Q yr -1 . Our analytically derived mass loss rates 
range from ^10~ u — 10~ 9 M yr -1 (when taking into 
account the error on the stellar flare frequency power 
law slope). For comparison, TTS wind mass loss rates 
have been derived from various line diagnostics: for TW 
Hva. lDupree et al.l (|2005l ) report M Q vn/0 = 2.3X10" 11 
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where 



M Q yr" 1 and M c m/</> = 1.3x10 
4> is the fractional stellar surface area from which the 
wind originates (e.g., a more collimated, polar wind will 
have cj> ^$0.3). These outflow rates are less than the Ha 
emission d erived accretion rate fo r TW Hya of4xl0~ 10 
M Q yr" 1 (jMuzerolle et all 120001 ). In cases of other ac- 
creting systems, wind mass loss rates inferred from ac- 
cretio n signatures in spect ral lines range from 10~ 9 to 
10" 7 (jHartigan et al.lfl995h . 
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4. ANGULAR MOMENTUM LOSS VIA CMES 

Having calculated mass loss rates from PMS stellar 
CMEs, we now focus on whether or how these episodic 
mass loss events could impact the angular momen- 
tum evolution of a TTS. We have determined order-of- 
magnitude CME mass-loss rates, calculated empirically 
as well as analytically, giving approximate lower and up- 
per limits of ^10~ 12 and ~10 -9 M Q yr _1 , respectively 
(Sections 13.2.11 and I3.2.2j) . In this section, we estimate 
magnetic lever arm lengths to calculate torques and as- 
sess their relative importance in spin evolution. 

4.1. What is the magnetic lever arm length? 

The torque on the star, due to a steady, magnetized 
wind can be written 

T w = -M w n*ri, (12) 

where M w is the mass loss rate in a steady wind, Q t is 
the angular rotation rate of the star, an d r\ is the aver- 
age "le ver arm" radius in the wind (e.g. JMatt fc Pudritzl 
2008a], and references therein). 

Computing the lever arm radius precisely requires 
knowledge of the three-dimensiona l flow structu re and 
magnetic held configuration (e.g., iMestell 11984) . For 
time- variable flows, such as CMEs, the situation is even 
more complicated. Thus, here we will only estimate the 
lever arm radius that is appropriate for an ensemble of 
CMEs and that is based upon our current understand- 
ing of magnetized stellar winds. Based on simulations of 
steady-state, ideal MHD winds from stars with dipolc 
magnetic field and rotating at 10% of breakup speed 
(jMatt fc Pudritzll2008at iMatt et al.|[20H found 



R* V McMEVcsc/ 

where is the equatorial field strength of the 

global/dipolar magnetic field at the stellar surface, i?„ is 
the radius of the star, and v csc is the escape speed from 
the stellar surface. Here, we have introduced a factor ?/, 
in order to account for the difference in the effective lever 
arm length for an ensemble of eruptive mass loss events 
(A/cme), as compared to a steady wind (M w ). 

There are several factors that modify the effective lever 
arm length in an eruptive outflow, compared to a steady 
one. The first is that the mass loss is not continuous, but 
happens in discrete bursts/ejections. For simplicity, we 
consider that the mass loss occurs with an average rate 
of M w = Mcme, but occuring via discrete bursts that 
happen for a fraction f t of the time (ft < 1). Thus, each 
burst has an instantaneous mass loss rate of / t _1 McME, 
and the factor of rj in equation [13] should include a fac- 
tor of ft, in order to take this into account. It is clear 
that a time-dependent wind is less efficient at extracting 
angular momentum than a steady wind with the same 
average mass loss rate. In the inferred distribution of 
CMEs (see Fig. [3J middle panel) , the mass loss is dom- 
inated by the lowest energy events, of which there are 
^10 per year. To estimate f tl we also need to specify 
the duration of each event — i.e., the duration over which 
we can think of the CME as a steady wind from the sur- 
face of the star. As an order of magnitude estimate, we 
consider the time for a CME traveling several hundred 



km/s to cross several stellar radii, which is of the order of 
one hour. For 10 events per year with a duration of one 
hour, ft ~ 10~ 3 . To calculate our upper limit to Mcme 
in Section f3.2.2[ we considered that the flare rate distri- 
bution extends to much lower energies than observed. In 
this case, the event frequencies of the lowest energy flares 
are much higher, implying a larger value for f t . 

The second factor modifying the effective lever arm 
length in an eruptive flow is that the CME events are not 
globally distributed in the corona, but take place over a 
fraction of the total solid angle of whole sky. More impor- 
tantly for the angular momentum flow, not all of the stel- 
lar magnetic flux participates in the azimuthal accelera- 
tion of the outflow. To take this into account, the factor 
r\ in equation ll3l should include a factor of where is 
the fraction of (unsigned) magnetic flux participating in 
the CME flow, compared to the total stellar surface mag- 
netic flux associated with the global/dipolc field. The un- 
signed flux participating in a CME should approximately 
equal the geometric area times the field strength of the 
underlying active region. Generally speaking, and al- 
most by definition, the magnetic field strengths of active 
regions are much stronger than the global field strength, 
while the areas are much smaller than the total stellar 
surface area. The product of these two approximates f^ , 
and it is not clear in general whether this will be very 
different from unity. For the purpose of this section, we 
assume ~ 0.1-1.0. 

Other factors that influence the effective lever arm 
length is the angle that the direction of the CME flow 
makes with respect to the rotation axis of the star and 
that the acceleration of the CMEs may be significantly 
different than the acceleration of the steady wind in the 
simulations of Matt et al. Given the much larger uncer- 
tainties discussed above, we neglect these effects for the 
present analysis and approximate r\ = ftfl- We consider 

here a range of 10~ 3 < rj < 1 but acknowledge that an 
even wider range may be possible. Note that the case of 
77 = 1 is equivalent to the case of a steady wind. 

To compute the lever arm length as a function of time, 
we use a model for the evolution of a 1 M Q star from 
iSiess et all (|2000f h which gives the stellar radius i?* (and 
thus surface escape speed v csc ) as a function of stellar age. 
The upper left panel of Figure Q] shows the evolution of 
i?*. To compute the lever arm length, we must specify 
the surface magnetic field strength, mass loss rate, and 
rj. For simplicity, we adopt a singular value of B* = 600 
G (corresponding to the equatorial field strength of th e 
dipole field measured on BP Tau; iDonati et al.l [20081 . 
and assume that it is constant in time. The range of 
possible values of for TTSs, as well as how this may 
evolve in time, is not well constrained observationally, 
which introduces an uncertainty in our calculated val- 
ues. The upper right panel of Figure [4] shows the Alfven 
radius as a function of time, for four combinations of the 
extreme values of mass loss rate and 77 discussed above. 
In the Figure, the vertical dash-triple-dotted lines mark 
the approximate age range of TTSs (for which this calcu- 
lation is valid), corresponding roughly to an age of a few 
Myr, plus or minus a few Myr. Since the mass loss rates 
we are considering correspond to TTSs, the calculation is 
only valid/relevant between these two vertical lines, but 
we show a wider range of ages for context/illustrative 
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purposes. It is clear that the range of rj and Mcme con- 
sidered corresponds to a range of a factor of about 13 in 
the possible Alfven radii at any given age. In the age 
range of TTSs, the Alfven radius could be in the range 
of 3-70 depending on the wind parameters. 

Observationally, the largest magnetic structures are 
seen to exist at distance scales compara ble to Alfven 
radii c alculated here: the magnetic loops of iFavata et al.l 
(2005) were inferred to be from 0.1-55 i?* in extent (as- 
suming these were indeed single loops), and cool promi- 
nences have been seen to extend up to a few stella r 
radii from TTS (|Massi et al.l [20(51 ISkellv et al.l l2008h . 
At later evolutionary stages, solar type stars are still 
seen to have extended, su spended material: the 50 Myr 
old (|Luhman et al.ll2005l ) solar analog AB Dor was ob- 
served to have cool, corotating clouds at distances from 
^9-20 R» from the rotati on axis, and were estimated 
to h ave masses >10 18 g (ICollier Cameron fc Robinson! 
119891) . iDunstone et~atl (|2006j) observed prominences 
on Speedy Mic, a 30 Myr old K3 dwarf, and esti- 
mated masses to be ~2x 10 17 g. In stellar outflows, the 
Alfven radii are typically comparable to or slightly larger 
than the size of the l argest closed magnetic loops (e.g., 
IMestel fc~S pruit 1987J). Thus, the measured sizes of post- 
flare loops and prominences may represent an approxi- 
mate lower limit on the Alfven radii for these systems. 

4.2. How important is the torque? 

The spin rate of stars will change in response to exter- 
nal torques and changes in the structure of the star (and 
subsequent change in the stellar moment of inertia). The 
conservation of angular momentum for an isolated star 
rotating as a solid body can be expressed as 



cK7» 
~dt 



T — fl 



dl* 

~dT' 



(14) 



where /* = k 2 M r R 2 is the moment of inertia of the star, 
Jl* is the angular spin rate of the star, T w is the wind 
torque (eq.[T2]); here, we're treating the CME outflow like 
a wind, so T w — Tqme), and k 2 is the normalized radius 
of gyration (defined by the radial mass distribution) . For 
a star in a tight binary system or a star that is still 
actively accreting from a disk, there will be additional 
torque terms in equation 1141 which we neglect here in 
order to isolate the effects of CMEs alone. 

To compute the relative s trength of t he ter ms in this 
equation, we use the 1 M Q Sie ss et ahl (|2000f ) model to 
specify the evolution of k 2 and In this and in most 
other PMS stellar evolution models, the star is fully or 
nearly fully convective during the first ~10 Myr of evo- 
lution. In convective regions, the mixing of material 
redistributes angular momentum on a timescale that is 
much shorter than evolutionary timescales. Thus, exter- 
nal torques effectively act upon the entire mass of the 
convection zone; for a fully convective PMS star, the 
external torques act on the entire star (as opposed to 
acting only on a thin shell that is rotationally decou- 
pled from the interior) . There have been observations of 
possible differential rotation in a few TTSs, with ampli- 
tudes of ~ 10% of the bulk/average rotation rate (e.g., 
iHerbst et all I2006H . but in the vast majority of TTSs 
studied no measureable differential rotation is observed. 
Thus here we follow convention and treat the star as a 



solid-body rotator, where represents the bulk rotation 
rate, and the error associated with this approximation is 
much smaller than for other unknowns (such as the value 
of rf). In order to specify f2*, we simply assume that 
the star is always rotating at a fraction of / = 10% of 
breakup speed, in order to approximately represent the 
"slow rotators" in observed TTS spin distributions (e.g., 
IHerbst et al. l l200l . This means that the rotation period 
is assumed to be 

Since both terms on the right hand side of equation H~4l 
are proportional to 57*, the relative value of these two 
terms is independent on the assumed value of 57*. 

The bottom left panel of Figure @] compares the two 
terms on the right hand side of equation Q3] as a func- 
tion of time. The broken lines show the value of — Tqme, 
where we have multiplied by —1 in order to compare 
the absolute value of the terms on a logarithmic scale. 
The different line styles and colors correspond to the as- 
sumptions and range of possible values of stellar wind 
parameters discussed in the previous section and labeled 
in the top right panel of the Figure. The stellar wind 
torque always acts to spin down the star (decrease O*). 
The solid line in the bottom left panel shows the last 
term on the right-hand side of equation [T3J This term 
has the dimensions of torque and can be thought of as 
a pseudo torque, which describes how the spin rate of 
the star changes due to changes in the moment of inertia 
(i.e., even if total angular momentum were conserved), 
and which acts to spin up the star (increase 57*). 

In order to assess how the terms in equation Q3] may 
affect the stellar spin rate, it is useful to compute a spin- 
up or spin-down timescale by dividing the total angular 
momentum by the torque. By rearranging equation 1141 
the spin-up/down time can be written 



-7* 



Tcme/I* - £l*/l*(dl*/dt) 



(16) 



When the CME-outflow torque dominates, the star will 
spin down on a timescale 



TCME 



ICME 



= k 2 
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(17) 
(18) 



By contrast, when the changes in stellar structure dom- 
inate, the star will spin up on a timescale 



^struct 



(dh/dt) 



(19) 



One advantage of looking at things in this way is that 
it is independent of the current spin rate of the star. In 
the bottom right panel of Figure HI we show T strU ct (solid 
line) and tcme (broken lines, corresponding to the vari- 
ous wind parameters indicated in the upper right panel), 
as a function of stellar age, for the 1 M Siess et al. 
model. At any given time, the shortest timescale is the 
dominant one, and if it is comparable to or shorter than 
the age of the star (shown as a dotted line) , it will have 
a noticible effect on the stellar spin rate. 
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It is clear from Figure 0] that the stellar contraction is 
expected to effectively spin the star up during the first 30 
Myr (in the absence of any external torques). The Figure 
shows that this pseudo-torque from stellar contraction 
completely dominates over the angular momentum loss 
associated with CMEs, for most of the range of parame- 
ters considered here. The only exception is that, for the 
highest mass outflow rates allowable by our analysis, and 
for T) ~ 1, the angular momentum loss by CMEs becomes 
comparable to the contraction pseudo-torque at an age 
of a few million years and dominates it thereafter. 

Many of the stars in the COUP sample are actively 
accreting from a surrounding accretion disk. It is there- 
fore interesting to consider the angular momentum ex- 
change between the star and disk , arising from this inter- 
action. iMatt fc Pudritz (2008b) demonstrated that the 
spin-down torque associated with the magnetic connec- 
tion between the star and disk will be less than the torque 
from a stellar outflow, in the case where the Alfven ra- 
dius is not too large and when the magnetic coupling to 
the disk is strong. Using their "preferred" values for the 
coupling, they demonstrated that a stellar outflow will 
carry away more angular momentum than the star-disk 
interaction, as long as ta < 84i?„». As demonstrated in 
the upper right panel of Figure IH ta < 70i?* for all of 
the outflow parameters considered here, so the spin-down 
torque due to a star-disk magnetic connection should be 
negligiblcQ- Accreting stars can also experience strong 
spin-up torques, mainly associated with the accretion of 
high specific angular momentum from nearly-Keplerian 
disks. The spin-up torque from accretion is given by 
M a y/GM*R t , where M a is the accretion rate, and R t is 
the radial location of the inner truncation o f disk (e.g., 
iGhosh fc Lalnbl[T97llM"att fc Pudritzll2005bfl . As an ex- 
ample, for the star considered here, accreting at a rate 
of lO^Afoyr" 1 , with a truncation radius of a few stellar 
radii, and at an age of 6 x 10 years (values appropriate 
for the early TTS phase), the resulting accretion torque 
is ~ 10 37 erg. This spin-up toruqe is comparable to the 
contraction pseudo-torque at the same age (lower left 
panel of Figured]). 

Thus, depending on the stellar age, accretion rate, etc., 
the processes acting to spin up the star may be dom- 
inated by either contraction or accretion. In order to 
counteract this spin-up torque and to explain the obser- 
vations of slow rotators at all TTS ages, the mass outflow 
rates must be high, the magnetic field strengths must be 
large, and/or additional spin-down torques (that are not 
considered here) must be present. The CME outflows 
inferred in the present work could be important near the 
end of the TTS phase, in the case of weakly- or non- 
accreting stars, and if the mass loss rate is near our upper 
limit (Mcme > 1O~ 1O M yr _1 ) and rj ~ 1. 

5. DISCUSSION 

In this work, we derive mass-loss rates for CMEs from 
solar type PMS stars. This is a first attempt to estimate 
PMS mass-loss via this mechanism and to calibrate it 
explicitly against the Sun and against observations of 
flaring magnetic loops on PMS stars. We find the to- 

1 The analysis of lMatt Pudritd l|20P8b) is valid only for r) = 1. 
It could be modified to consider 7] ^ 1, but the overall results are 
not affected. 



tal CME mass- loss rate to be in the range of 10~ 12 to 
10 -9 M Q yr" 1 . This is relatively modest by compari- 
son with mass loss rates observed for TTS winds. How- 
ever, the magnetic lever arms associated with the largest 
PMS CMEs, as inferred from the COUP sample stars 
and scaled from the Sun, are large — up to tens of stellar 
radii — so the resulting torques can be significant. Assum- 
ing the CME mass loss rate for a typical ~1 M© TTS 
is < 10~ 10 M Q yr" 1 , the associated spin-down torque is 
likely too small to counteract the spin-up effects of con- 
traction and accretion in the T Tauri phase. However, 
the spin-down torque from CMEs could be important af- 
ter an age of a few Myr if the mass loss rate is > 10~ 10 
M Q yr _1 , the accretion rate is low, and the angular 
momentum loss is efficient. Our estimates of the an- 
gular momentum loss efficiency (i.e., the magnetic lever 
arm radii; see Section |4]) are based upon axisymmetric, 
steady-state wind calculations, and 3-dimensional simu- 
lations are needed to improve these estimates. 

Our stellar CME model is derived from relating solar 
activity to stellar; in one sense, we immediately under- 
estimate mass losses via CMEs by only selecting solar 
CMEs that occurred with flares. We find, however, that 
much of the mass loss comes from the CMEs associated 
with flares. The most massive CMEs are associated with 
the most powerful flares, and the most e nergetic flares 
are t he most often associated with CMEs (| Aarnio et al.l 
I2TJT1 . Thus, using stellar flare activity to determine a 
stellar CME rate likely represents a lower bound on the 
stellar CME activity. 

One of the most pivotal parts of this calculation is the 
extrapolation of the CME mass/flare energy relationship. 
The adopted direct extrapolation, while crossing many 
orders of magnitude, is indeed physically motivated, and 
introduces the fewest new assumptions. We are able to, 
using the COUP "superfiare" observations, justify the 
extrapolation up to TTS flare energies; our premise is 
that these loops represent proxies for stellar CME masses 
as the loops exist within the density and height regimes 
from which CMEs are generally launched on the Sun. 
The consistency with the solar CME mass/flare energy 
relationship is particularly compelling given that it spans 
six dex in flare energy. 

In ad opting the stellar X-ray flare fr equency relation- 
ship of (|Albacete Colombo et al.ll2007l ). we have gener- 
ated an "ensemble mass loss rate" ; in doing this, we have 
neglegted to separate out any mass-dependent charac- 
teristics of the distribution. There would be, however, 
significant differences in the flare activity for stars of the 
same mass with different rotation rates. Given the mixed 
sample of classical and weak-line TTS in the COUP, 
there are also self-absorption effects (by the accretion 
columns, iGregorv et al.l 120071 or a circumstellar disk 
inclined to our line of sight) which could introduce bias 
into a flare frequency distribution. It was our intent in 
using the COUP sample to diminish all of these obser- 
vational effects via the large number of stars observed in 
various activity and evolutionary phases, as well as in- 
clinations and rotation rates. Future work could include 
analyzing subsets of the COUP dataset to see whether 
there is an X-ray luminosity dependence on the flare rate 
(and thus, a difference in CME mass loss rates as well). 

As noted in Section 15 .2. 21 one particular difficulty with 
this analysis was the sensitivity of our analytical solu- 
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tions to the value of the power law slope. Within the 
errors on these quantities, we find our mass loss rates 
derived analytically could span two orders of magnitude. 
The empirically derived mass loss rates, while derived 
from incomplete samples, represent robust lower limits. 
In the solar case, our analytical mass loss rate did pre- 
dict more mass loss than was observed in the LAS CO 
database, but as discussed in Section [2 the LASCO 
database itself is missing many CME masses, so the to- 
tal mass loss rate is undoubtedly higher. Despite this 
apparent overprediction, we still derive a mass loss rate 
below the steady solar wind mass loss rate, which is as 
anticipated. 

Observing stellar CMEs would be ideal for address- 
ing the most uncertain parameters in this work, but 
presently, understanding the observational signatures 
and being able to detect a stellar CME remain outstand- 
ing probl ems. Efforts hav e been made with EUV spec- 
troscopy; iLeitzinger et al.l (|2011l) point out a lack of si- 
multaneous solar spectra during flare/CME events, cre- 
ating great difficulty for the interpretation of stellar data. 
Additionally, managing to catch a stellar CME with the 
right set of parameters for observation (e.g., density and 
projection) requires sufficient observing time to increase 
the probabilit y of detection. F o r a st eady, hot (1 MK), 
coronal wind, iMatt fc Pudritzl (|2007l ) calculate a mass 
loss rate of lxlCP 9 M Q yr -1 would produce strong EUV 
and X-ray emission, much higher than is observed. Our 
highest mass loss rate estimate, ~10 -9 M Q yr -1 , is prob- 
ably unlikely given the strong observational signature 
that would attend such a mass loss rate. Interestingly, 
however, the authors also found that decreasing the mass 
loss rate by a factor of 10 decreased the excess emis- 
sion by a factor of 100. Given the "bursty" nature of 
our CMEs, and the likely high cadence of these bursts, 
and potentially cooler temperatures (solar CMEs are ob- 
served to rapidly expand as they travel away from the 
Sun) , a significant level of mass loss could easily escape 
detection. 

The COU P "superflaring" tar gets were generally weak 
lined TTS (|Aarnio et al.l l2010f ); it is likely that accre- 
tors could show an enhanced activity rate and thus 
more frequent CMEs, a s well as accretion-driven winds 
(jMatt fc Pudritdl200 8a) which could further deplete an- 
gular momentum. The presence of a close-in disk in ac- 
creting systems might mitigate the enhanced CME rate 
with reduced lever a rm lengths, however modeling work 
(|Orlando et al.|[20Tll) has shown star-disk flaring can dis- 
rupt circumstellar disk material, causing an MM insta- 
bility that then generates an accretion flow. Thus, in a 
cyclical fashion, large scale flaring can lead to accretion, 
the accretion then powering winds, outflows, and more 
reconnecti on, result i ng in further mass / angular momen- 
tum loss. iCranmeil (|2009f ) illustrate one mechanism by 
which accretion could power a wind: infalling mate- 
rial impacting the stellar surface would generate Alfven 
waves, the propagation of which could trigger reconnec- 
tion, accretion thus powering activity. 

6. SUMMARY AND CONCLUSIONS 

This represents a preliminary effort at calculating mass 
and angular momentum losses for T Tauri stars via 
scaled-up, solar-analog coronal mass ejections. Our cal- 
culations are based on well-observed flares among a large 



sample of solar-type pre-main-sequence stars, and our 
analysis methods have been tested for the most well- 
understood case, the Sun. Beginning with an established 
relationship betwee n solar flare energies and CME masses 
(| Aarnio et al.ll2011f ). we use the observed flare frequency 
distribution to infer a CME mass distribution. We are 
able to infer empirically and analytically solar CME mass 
loss rates that are consistent with observations, and also 
consistent with expectations for a mass loss rate once 
observational biases are carefully accounted for. 

Our lower limit on the mass shed via CMEs in a generic 
TTS case is 6.2xl0~ 13 M yr" 1 . To obtain an up- 
per limit on the stellar CME mass loss rate, we calcu- 
late analytic solutions to account for observational bi- 
ases (most importantly, the stellar X-ray flare detection 
limit); our upper limit to the mass loss rate via CMEs is 
M C ME<2xlO^ 9 M yr- 1 . 

Finally, we assess the resulting torque against the 
star's rotation provided by the CME mass loss. We find 
that within the ranges of mass loss rates and effective 
lever arm lengths (reflecting the difference in steady ver- 
sus "bursty" mass loss) being considered, the resulting 
Alfven radii would span a large range, 3-70 R*, con- 
sistent with previous inferences of the physical sizes of 
large-scale magnetic structures. We find that near our 
upper limit mass loss case (>10 -10 M Q yr _1 ), this mech- 
anism could be effective for slowing stellar rotation on a 
timescale comparable to TTS lifetimes. 

We acknowledge NSF grant AST-0808072 (K. Stassun, 
PI). 
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Figure 1. Solar relationship between flare energy and CME mass, transformed from the flux/mass relationship of Aarnio et al,l II201H ). 
The line is a fit to the solar data over the energy range of observed solar flares, 10 26 -10 30 erg, and then extrapolated to the regime of T 
Tauri star flare energies (represented by vertical l ines). The single po int at upper right shows the mean CME mass and flare energy for the 
32 T Tauri star mega-flares observed by COUP (Favata ct al. 2005), and the error bars represent the standard deviations of those mean 
values. The extrapolated solar relationship fits this point well. (Inset:) The energy distribution of flare rates for T Tauri stars observed by 
COUP (Albaccte Colombo ct al. 2007). The vertical dotted line at 10 35 6 erg represents the energy above which the observed flare sample 
is complete. 
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Figure 2. Observed distributions of solar flare energy and of m easured solar CME masses, compared to CME masses inferred from the 
observed flares using the methodology described in Section |3. 1.11 In the left column, flare rates as a function of energy from the GOES 
database from 1996-2006 are shown (histogram, upper panel; differential dN/dE distribution, lower panel). The middle column shows 
CME rates as a function of masses from the LASCO database during the same period (histogram, upper panel; differential distribution 
dN/dM, lower panel). In all four of the observed distributions (left and middle columns), black curves denote the full distribution, and red 
dot-dashed curves show the distributions of properties for only those flares and CMEs which are associate d. The CME mass distribution 
includes all CMEs in the LASCO database with reported masses plus the Halo CMEs assigned masses bv lAarmo et al.1 (2011 ) ■ Panels m 
the right column show our inferred CME rate as a function of mass, inferred from the observed flare distributions described in Section [XT] 
The black distributions are the inferred distributions, while the gray are the observed distributions (as seen in black in the center column), 
for comparison. In bottom panels, vertical dotted lines represent the energy or mass limits above which the distributions are approximately 
complete. 
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Figur e 3. Flare and CME d istributions for PMS stellar sample from the COUP. The left panel is the flare frequency distribution of 
Albacctc Colombo ct al. (2007), expressed in terms of number of flares of a given energy per time per energy bin (dN/dE). In the middle 
panel, we show our inferred TTS CME mass distribution: we have taken the product of dN/dE from the left panel and dE/dM derived 
here in Section[2]to arrive at dN/dM. Then, in the right panel, we show the stellar CME mass distribution were it observed over the same 
time span as the LASCO CMEs (compare to Fig. [2] upper middle panel). 
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Figur e 4. Uppe r left panel: Evolution of the photospheric radius of a one solar mass star, as a function of age, from the model tracks of 
ISiess et al.l 112000) . Upper right panel: Effective magnetic lever arm (Alfven) radius in CME outflows from the same star, with a magnetic 
field strength of 600 G and four different combinations of mass loss rate (Mcme) and efficiency factor as indicated on the plot (Mcme 
is given in units of Mq yr - 1 ). The vertical dash-triple-dotted lines show the approximate age range of the T Tauri phase. Lower left 
panel: The solid line shows the "pseudo- torque," which shows the contribution to the evolution of stellar rotation rate that is due to the 
contraction and changes of internal structure of the star(see text). The colored lines show the CME-outflow torque, which acts to spin 
down the star, corresponding to different outflow parameters (line styles correspond to the upper right panel). Lower right panel: The 
solid line shows the spin-up timescale, due only to the contraction of changes of internal structure of the star, as a function of stellar age. 
The colored lines show the spin-down timescale due to CME outflows, corresponding to different outflow parameters (corresponding to the 
upper right panel). The dashed line corresponds to a spin-timescale equal to the age of the star. 



